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Abstract 


We consider a static self-gravitating charged perfect fluid system 
in the Einstein-Maxwell theory. Assume Maxwell’s equation and 
the Einstein constraint equation are satisfied and the temperature 
of the fluid obeys Tolman’s law. Then, we prove that the extrema 
of total entropy implies other components of Einstein’s equation for 
any variations of metric and electrical potential with fixed boundary 
values. Conversely, if Einstein’s equation and Maxwell’s equations 
hold, the total entropy achieves an extremum. Our work suggests 
that the maximum entropy principle is consistent with Einstein’s 
equation when an electrostatic held is taken into account. 

PACS number(s): 04.20.Cv, 04.20.Fy, 04.40.Nr 

1 Introduction 

It is well known that black holes can radiate and satisfy thermodynamical 
laws. This discovery establishes the connection between gravity and ther¬ 
modynamics [1-6]. Despite the great success of black hole thermodynamics, 
there are still unresolved issues, for instance, the origin of black hole en¬ 
tropy. In contrast to black hole systems, local thermodynamic quantities of 
a perfect fluid in curved spacetimes, e.g., energy density p, entropy density 
s, and local temperature T, are well defined. The presence of gravity only 
affects the distribution of those local quantities. In general there are two 
methods to determine the distribution of matter. One way is solving Ein¬ 
stein’s equation. The other way is using the entropy principle to determine 
the distribution of matter. 

Since entropy plays no role in Einstein’s equation, it is unclear whether 
the two methods are consistent. Even before the establishment of black 
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hole thermodynamics, Cocke [7] pointed out that the extrema of entropy 
should yield the equation of hydrostatic equilibrium which is derived from 
Einstein’s equation. After that, Sorkin, Wald, and Zhang [S] showed rigor¬ 
ously that the Tolman-Oppenheimer-Volkoff equation of hydrostatic equi¬ 
librium can be derived from the extremum of total entropy and the Einstein 
constraint equation. Gao [5] extended their proof from radiation to a gen¬ 
eral perfect fluid, including uncharged fluid and uniformly charged fluid. 
This issue has been further explored in the past few years [10-18]. 

Recently, we m proved the entropy principle for a self-gravitating fluid 
in static spacetimes without any symmetry in the spacelike hypersurface. 
So far, the matter field considered is a perfect fluid. It is interesting to 
know whether or how the entropy principle works in the presence of an 
electromagnetic field. In this paper, we extend the two theorems in Ref. 
[19] to a uniformly charged fluid. The extension is not straightforward at 
all. For an uncharged perfect fluid, we have shown [12] that the variation 
of total entropy, SS, is proportional to the variation of the spatial metric, 
hab- However, for a uniformly charged fluid, 6S appears to depend on 6hab, 
dx, and 6Aa, where x is the redshift factor and A“ is the vector potential of 
the electrostatic field. However, we manage to show that the components 
of Sx and Shab vanish identically and the vanishing of the components of 
5hab just gives the spatial components of Einstein’s equation for a charged 
perfect fluid. Our work suggests that the entropy principle is consistent not 
only with the gravitational field but also with the electromagnetic field. 

2 Properties of charged perfect fluid in static 
spacetimes 

We consider a general perfect fluid as discussed in Ref. |9]. The entropy 
density s is taken to be a function of the energy density p and particle 
number density n, i.e., s = s{p,n). From the first law of thermodynamics, 
one can derive the integrated form of the Gibbs-Duhem relation, 

s = ^{p+p-f^n), ( 1 ) 

where p and p, represent the pressure and the chemical potential, respec¬ 
tively. All the quantities are measured by static observers with 4-velocity 
u°‘. These observers are orthogonal to the hypersurface E. Therefore, the 
induced metric on S is given by 


hab — Qab UaUb ■ 

The stress-energy tensor Tab for a perfect fluid takes the form 

Tab = pUaUb -I- phab = {p + p)UaUb + PQab ■ 


( 2 ) 

(3) 
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Since the fluid is charged, we should also consider the stress energy 
tensor of the electromagnetic field, 

{T^^)ab = ^ {PacFb^ - \gabFdeF'^<^'^ , ( 4 ) 

where Fab — Vq^{,—V hAa, and Aa is the vector potential. The electrostatic 
potential $ is given by 

<^> = -AaC = -AaU‘^X, ( 5 ) 

where is the Killing vector and x is the redshift factor. 

The Maxwell’s equation with source reads 

= 47rj“ = 47rpeM“ , (6) 

where and Pf, represent the 4-current density of the electric charge and 
the charge density, respectively. Then 

^'^Fab = ^\VaAb-VbAa) 

= ^ai^'^Ab) - ifeVaC"' - C'Vfeia 

= -Va4>, (7) 

where we have used the Killing equation and C^Aa = 0 in the last step. 
From VjaFfjc] = 0 we can calculate 

y^(j.EM)ah ^ 

= -peU‘=F\ 

X 

= ( 8 ) 

X 

where Maxwell’s equation S/[aFbc] = 0 has been used in the first step and 
Eq. 0 has been used in the last step. 

We shall assume that Tolman’s law holds, which states that the local 
temperature T of the fluid satishes 


Tx = To , (9) 

where x is the redshift factor for static observers and Tq is a constant. 
Without loss of generality, we take Tq = 1. This law establishes the rela¬ 
tionship between the fluid temperature and the metric components. 

It is then straightforward to show, from the conservation law Vq [T“^ -|- 
^rpEM'^abj _ Q stationary conditions, that 

0 = V>-f (p-f -h — V“4>, (10) 

X 
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where A°‘ is the 4-acceleration of the observer. For stationary observers, 


Aa = ^aXlX-, 

(11) 

and thus 


VaP=-(p+p)VaX/x-^V“<I>. 

X 

(12) 

On the other hand, the local first law can be expressed in 

the form [9] 

dp = sdT ndfjL. 

(13) 

Using Eqs. (O and (IT|), we find 


VaP = + 


X 

(14) 

Comparison with Eq. (IT^ gives 


- —= AiVaX + XVaM = Va(/ix) ■ 

n 

(15) 

If we assume that all particles possess the same charge q, i.e.. 


Pe = nq, 

(16) 

then Eq. (fT^ leads to 


^X H- 9$ = c, 

(17) 

or 


- -h 9$ = c, 

(18) 


where c is a constant. 


3 Two theorems 

The distribution of a charged perfect fluid in static spacetimes can be deter¬ 
mined in two ways. First, Einstein’s equation and Maxwell’s equation to¬ 
gether can totally determine the distribution. Second, under certain bound¬ 
ary conditions, the total entropy of the fluid should take an extremum. In 
this section, we prove two theorems indicating the equivalence of the two 
methods. 

Theorem 1: Consider a uniformly charged perfect fluid in a static 
spacetime {M,gab) and E as a three-dimensional hypersurface denoting a 
moment of the static observers. Let C be a region on E with a boundary 
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C. Let hab, and x be the induced metric on E, electrostatic potential, 
and redshift factor, respectively. Assume that the temperature of the fluid 
obeys Tolman’s law and the Einstein constraint equation and Maxwell’s 
equation are satisfied in C. Then, the other components of Einstein’s equa¬ 
tion are implied by the extrema of the total fluid entropy for all variations 
of data in C where hat, Xi and their first derivatives are fixed on C. 

Proof - The total entropy S is an integral of the entropy density s over 
the region C on E, 


S = [ Vhs{p,n), (19) 

Jc 

where h is the determinant of hab- Thus, the variation of the total entropy 
is written in the form 


6S= sS'/h + Vh6s. ( 20 ) 

JC 

Applying the local first law of thermodynamics, 

Tds = dp — fidn , (21) 

we find 



Note that p/T is constant for an uncharged fluid and can be moved out of 
the integral [19]. But due to the electrostatic potential, Eq. m shows that 
p/T is no longer a constant. We shall deal with the 5n term by employing 
Maxwell’s equation. 

Together with Eqs. CD and CHD, we have 

5S = [ —{p + p—pn)6Vh—Vh{—q^ + c)6n + Vh—Sp. (23) 

Jc ^ ^ 

Denote 

6S= f SL, (24) 

Jc 

where 

5L = —{p + p— pn)6Vh — Vh{—q^ + c)Sn + Vh—Sp . (25) 

Our purpose is to derive the space components of Einstein’s equation 
from 6L = 0 and the constraint Einstein equation. First, we need to express 
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6L as variations of basic variables hab, X s-nd $. The 5h term in Eq. (1^51) 
can be easily written in the desired form by the relation 

6y/h=^Vhh‘^^Shab. (26) 

The Sn term in Eq. (051) is calculated by using Maxwell’s equation (see 
Appendix B). Now we shall focus on calculating the 6p term. 

Note that the extrinsic curvature of E defined by 

Bab = Khi^dUc (27) 

vanishes [19] in static spacetimes and 

VbUa = -AaUb. (28) 

By the result of Ref. m, the Ricci tensor and scalar curvature 
of E are given by 

^l^b = ^ab + RaebU'^Ul + RfbU^Ua + RakU^Ub + UaUbRfkU^ (29) 

and 

i?(3) = . (30) 


To calculate 5p, we start with the Einstein constraint equation 

Gabu’^v!' = 87rT*f , (31) 

where 

= + (32) 

Together with Eqs. and m, we find 

^ ^ (33) 

Denote the last term of Eq. (1251) by 6Lp. By substituting Eq. (p3)) into 
Eq. (I25|), we have 


SLa = 


1 


-VhSR^^'> - Vh 

IGttT AttT 


SiFacFb^U<^U^) + -SiFabF^^^) 


(34) 


By the standard calculation [TH], the first term on the right-hand side of 
Eq. (1341) can be written in the desired form 


1 

IGttT 


VhSR^^'> 


— (-^-R(3)ab J_Mf)Shab , 

T ^ Wir IGtt ^ ^ ’ 


(35) 
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where 


+ D^A°' - 


(36) 


and Da is the derivative operator associated with hab- Denote the second 
term on the right-hand side of Eq. (IMl) by SLp 2 , i.e., 


6L 


p2 




AttT 


5{FacFb^U^U^) + ^5{FabF^^) 


(37) 


The calculation of Eq. (l37)) is given in Appendix A. From Eqs. (|60l) and 
|), we have 


^Lp2 — 


2Da{u^Fb^)5<^ + U^Fa\D‘^'i>)5hab + 


'/h 

At: 


Vh 

At: 
y/h 
At: 


X 


2Da{u^Fb^) - AaF-^Ub + Vai^) — 

T X. 


i5<i) 


u‘^Fc\D^<A?) + 


^hab 


VhDaA>D^A> 


Attx^ 


5xm 


The substitution of Eqs. (|35|l and (l38l) into Eq. (IMl) yields 


SL, = Ai 

H rji 


-Li?(3)a& jjb _ J^aby 

IGtt IGtt 


^^ab 


+ 


y/h 

At: 

yfh 

At: 


, Ub 


2Da{v!>Fb‘^) - AaF^\b + Va( V) 

1 




^hab 


T ' X 


At:x 


5(^ 

Cfb 

2 ^X-(39) 


The substitution of Eq. (15^ and Eqs. (1751) and (1751) in Appendix B 
into Eq. (Eg yields 


5L = 


^{p + p — pn)6y/h — y/h{—qA> + c)Sn + y/h^Sp 

SL^ -|- SL^ F SLfi , (40) 


where 


6L^ = 


- — y/hD\u‘^Fcb)6<A> 
47r 

y/h 


+ ' 


47r 


E'*", Ub 


2Daiu‘’Fb^) - AaF^^^Ub + Va( —) 


T ' x\ 


5^, (41) 
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(42) 


5L^ = 


Vh 2D‘=$L»c4> 


4tt x^ 


-<5x, 


SLh = + P ^'^ h'^Hhab - ^Vh{D‘^^)x-^D'’^5hab 


Vhh^'^Dd 


T 2 
1 

Stt 

\/h 


47r 

($ - 

<7 


^hab 


-Li?(3)ab jjbj^a _ 

IGtt IGtt 


+ ■ 


'/h 

47r 




6h. 


■ab • 


(43) 


It is obvious that 5L^ = 0, which shows that the variation of the redshift 
factor X has no contribution to 6S. Now we show that 6L^ also vanishes. 
We calculate 


Da{u^Fb^) - AaF‘^‘’ub + Va 


6L^ = -^'/hD\u‘^Fab)6<^ + ^ 

47r 47r 

y/h 

Att 

y/h 

Att 

yfh 

Fk 


,b T?. c^ 


yab. 


1 cr 

Ub 

V T y 

X. 


i5<i) 


1 

S’ 

1 cr 

Ub 

V T y 

X. 


(5$ 


hWaiu'^Fb^) - AaF^^’ub + ^“ ( ^ 7 ) - (7 


yab 


SF 


Tpab 

,b jp c\ I ^,a^, V7 A Z 7 'Q, 6 ^,^ I V7 {T?o.h^, \ ^ 


Vciu'^Fb^) + u“UeVa(u'’F6") - AaF‘^‘>Ub + Va(F“'>Ub) - 


T 


Here, we have used Tx = 1. The last term of Eq. (HI) vanishes because 

y/h (VaUb)x - UbVaX 


4tt T 
V^™b/ 


X 


= -p-F°-^{-AbUa- AaUb) 

47r 

= =0. 


(45) 


Note that 


(5$. 
(44) 


Valtc = -UaAc , 


(46) 




















and hence 


^ ['^ciu^Fb^) + U^UcVaiu^Fb^) - AaF^'^Ub + Va(F“'’ufa)] 6^ 

= ^ - A^F-’^Ub] 

= ^ [uVb'^^aU^ - AaF'^’^Ub] ( 5 $ 

Att '- ■' 

= ^ [-uVb^UaA^ - AaF^^’^Ub] d<^ 

Att 

= 0. (47) 

This result reveals that the variation of the electrostatic potential $ has 
no contribution to SS. So we have 


SL = SLb 


^ P + P _ ±^(^n-,^)x-^D^<i>Shab 

1 Z 47r 

+ ^Vhh‘^^Dd 

8tt 


$- ] X 

q 


^hab 


\/h 

y/h 

Att 


-Li?(3)a& + Jjbj^a _ /^aby^^c) 

IGtt IGtt 


u'=Fc'’(£>“$) + ^F^cF^^ 


Sha 


(48) 


This shows explicitly that 6S is determined by the variation of hab only. 
Since SS = 0 by the assumption of Theorem l,we have 


'^P±P^A^,ab 
T 2 


y/h 

47r 


OTT 




J_^(3)ab ^b^a 

IGtt IGtt^ 


/l“^VcA“) 


+ ■ 


y/Ti 

47r 


^c^^b(^a^) _pa 


,F' 


be 


= 0 . 




(49) 


By substituting Eq. (I55|) and Eqs. (1751) and (1751) in Appendix B into Eq. 
((55D . and letting c' = —cjq, we have 

+2r($ + c')h‘^^UcVdF‘^’^ - 2Th“''L»e[($ + c')udF‘‘‘^] 
+AT{X~^D^‘A>)D^<A> + _ (^a^b pbj^a _ 

-dTu'^Fc- 2h^^h^'^FceFd‘'. (50) 
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From Ref. [TO] we already know that 


+ D^’A’^ - = 0 . (51) 


Substituting Eq. (ED into Eq. (1501) . we have 
87rp/i“^ = h‘^^h’’‘^Rcd - 

+P“'’ + Pf , (52) 

where 

= 2FdcFe‘'u^u’^h^^ + 2r($ + c')h‘^^Uc'VdF‘^’^ 

-2r/i“^Pe[($ + c')udF‘‘‘^], (53) 

pf = 4r(x-ip''$)p“$ - 4ru=Pc'’£'“^> • (54) 


Now, we show that Pf ^ and Pf ^ vanish. We first calculate 


Dc{udF'^-) 


K^^eiudF’^^) 

U,U^Ve{udF'^^)+V,{udF‘^^) 
UcU^F‘^^VeUd + F'^^V.Ud + UdVcF^^^ 
-U,U^F<^‘^UeAd - F'^'^U.Ad + UdVeP^" 


(55) 


With the help of Eqs. © and Eq. (1531) can be written as 

P“'> = /i“^[2r2v'=$Vc$ + 2P(4> + c')ucVdP''^ 

-2P(4> + c')DciudF‘^‘^) - 2TudP‘^“Pc(4' + c')] 

= 2/i“'’[T2V'=4>Vc$ + r(4> + c')McVdP“‘' 

-r(4> + c)ud^cF'^‘' - T^V“$Vc4>] 

= 0. (56) 

Applying Eqs. 0 and ([9]) again, we find immediately 

Pf = 0 . (57) 

Therefore, Eq. (l58l) just gives the projection of Einstein’s equation on S 

Snph-^ = - 2(P,eP/ - \gcdFefF-f)h<^-h^'^ . (58) 

This completes the proof of Theorem 1. 

In the above proof, we used the Einstein constraint Eq. (1551) to derive 
Eq. (l25l) . Then, by applying 5S = 0, we obtained the spatial components of 
Einstein’s equation. It is not difficult to check that the proof is reversible; 
i.e., from the projected Einstein equation (64), one can show (5P = 0 in Eq. 
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(Hil), which makes the total entropy an extremum. Thus, we arrive at the 
following theorem. 

Theorem 2: Consider a perfect fluid with charge in a static spacetime 
{M,gab) and S as a three-dimensional hypersurface denoting a moment of 
the static observers. Let C be a region on E with a boundary C, hab, 
and X be the induced metric, potential, and redshift on E. Assume that the 
temperature of the fluid obeys Tolman’s law and both Einstein’s equation 
and Maxwell’s equation are satisfied in C. Then, the fluid is distributed 
such that its total entropy in C is an extremum for all variations where 
hab, ‘h, X, and their first derivatives are fixed on C. 


4 Conclusions 

We have rigorously proven the equivalence of the extrema of the entropy 
and Einstein’s equation under a few natural and necessary conditions. Dif¬ 
ferent from the proof for the uncharged case, we have to consider the vari¬ 
ations of X and A°‘. The treatment of the Sn term also totally differs from 
the uncharged case because g/T is no longer a constant. The significant im¬ 
provement from previous works is that we extended the maximum entropy 
principle to the Einstein-Maxwell theory. Our work suggests a clear con¬ 
nection between Einstein’s equation and the thermodynamics of a charged 
perfect fluid in static spacetimes. 
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A Calculation of Eq. (I37l) 


In this Appendix, we will show the detailed calculation oi — -^;S{FacFb‘^u°'u^) 
and — ^j, 5{FabF°'^) in Eq. (l37l) . First, we calculate 


y/h 
y/h 
y/h 


5{F^,Fb^u<^u^) 

[u^Fb^diFacU^) + U^Fac5{Fb^u'^)] 
[u^Fb^5{x~^D,<^) + 


47rr 

y/h 

'Iff 

y/h 

Fk 


[x-^v^Fb^5{D,<^) + />Fb^D,<^5x-^ + x~^u'^Fac5{h<^‘^Dd<^) + 


X-^/’Fb^DM - 




5x + x w“FacIi Fd(5$ + X u°-Fac{Dd‘^)5h 


dc 


u^Fb^DJ‘^> - + u^Fach^^^DdS^ + u^FaciDd^)6h'‘^ 

X 


(59) 


where Eqs. © and m have been used. Using integration by parts and 
dropping the boundary terms, we have 


y/h 

~A/f 

y/h 
Fn 
y/h 
47r 


5{FacFb^U^v!^) 

D/v!>Fb^)5F + + Dd{,u‘^Fach'^/6F + u<^F/{D’^^)5hdc 


X 


2D//‘Fb/5<F + ‘^E^^Sx + u‘^F/{D<^<F)5hab 
X 


(60) 


Now we turn to —-^^5{FabF°‘^). 

- ^^6{FabF'^’’) 

IGttT ^ ^ 

= -^^{F‘^^SF,b + FabSF‘^f>) 

= [F'^^SFab + F„fc5(F'=V"5“)] 

= [F‘^^SFab + Fabg^^g^‘^SF^‘‘ + Fabg'^'^5{F<^'^g'>'^) + Fabg^'^5{F^'^g'^/] 

= --^(2rtF„6 + 2F,,F6^(5g“*'). (61) 

iOTTi 
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Because SFab = 2V[a(5^b] and g°‘^ = h°-^ — u°-u^, we have 
'/h 


IGttT 
y/h 
IQttT - 
y/h 
IGttT - 
yfh 
IGttT - 
yfh 


6{FabF^f 

4F“'’V[a<5A] + 2FacFb^5h‘^'‘ - 2FacFb‘^5{u‘^v!’) 

4F“'’VaM6 + 2FacFb^Sh‘^^ - 2FacFb^u‘^5f - 2FacFb^fSu° 


AF^'^WaSAb + 2F„,Ff,"<5h“'> - 2FacFb^uTSx~' - ^F^cFb^F 


_ cciab 


c^a^b —1 


c^.l 


AD 

AF'^'^VJAb + 2FacFb^Sh^^ + 


167rT 

Since Aa = —^{dt)a, we hav^il 


SAa = -idt)aSA> = . 

X 

Using integration by parts for the first term in Eq. dnn), we have 

^ F^^X/aSAb 


AnT 

yfh 

'aft 


Faby / 


X 


- -^V 

— . V a 

47r 

y/h. 


ub 

~T 


5$ 


yfh^ 

A-k 


F^^\ Ub 


T J X 


-SA> 


yfh 


T'ah 


= -FvXF‘W4.) + ^VjU-)Hij<I.^ 


47r ^ ^ 47r V ^ / X 

Note that for any tangent to S, i.e., v°‘Ua = 0. Thus, 

DaV<^ = ha^hd‘^V,v‘^ 

= yaV‘^-v‘^Aa, 


or 


Therefore, 


= DaV^ + V^Aa . 


Att 

= -^DfF'^’^UbS^) - ^AaF‘^^Ub5^ 

Att 47r 

= -^AaF-^UbSA>, 

Att 


^ The coordinates are fixed for variations. So 6{dt)a = 0 


(63) 


(64) 

(65) 

( 66 ) 


(67) 


'r<5x-' 

(62) 
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where we have dropped the boundary term in the last step. Then Eq. (IMl) 
becomes 


y/h 




47r 47r 


ah 


F' 

~ 


X 


(5$. 


( 68 ) 


The substitution of this result into Eq. (l62|) gives 


y/H 

IGttT 


5{FabF-^) 


y/h 

IGttT 


4D 

4F“VaMb + 2F„eF,‘=h“'’ + ^4- 6x 

X 


y/h 


= -^AaF^\b5^ + 


y/h 


Fn 

y/h 

'%/f 


47r 


F' 


ah 


T 


Ub 


, D 

F\F^^6hab- " ^ 6x. 

Anx 


(69) 


B 


Calculation of the 6n term in Eq. 


m 


In this Appendix, we will show the detailed calculation of —y/h{—qF + c)Sn. 
Define 

SSn = — I y/h{—qF + c)6n = j SLn ■ (70) 

Jc Jc 

To express Sn as a combination of (5%, <54), and Shab, we need to employ 
Maxwell’s equation. Suppose that all particles possess the same charge q, 


Pe = qn. 


From Maxwell’s equation Eq. (jG)), we have 

n = -^UaVb(F“''). 
Airq 

Therefore, with the help of Eq. (IS51) 


n 


-—UaX/bF‘^^ 

47rg 

-^VbiUaF^^/ + 

47rg 47rg 

-^Db{UaF-/ - ^UaF-^Ab 
47rg 47rg 

-^^Dbix-^D^F). 


^F‘^\bAa 

47rg 


(71) 


(72) 


(73) 
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Since 


T 


= -q^ + c, 


where c is constant, we have 


6Ln = —Vh^—q^ + c)6n 

y/h 


Airq 


{-q^ + c)5[Di,{x 


Note 


( 74 ) 


( 75 ) 


Db{x~^D^A>){X) = Dbmx~^h!>‘'D,A^) + C\d{X)x-^D'^^, ( 76 ) 
and then 

5Ln = ^^{-qA^ + c)5[D,{x-^D^A^)] 

= ^(-<7$ + c)D,{D'’A^5x-^) + ^(-<7$ + c)D,{x-^DH^) 

Airq Airq 

+ ^(-g$ + c)Db{x~^DMh^l + ^(-g$ + c){x-^D‘^^)5C\d 

Airq Airq 

= ^ 74 ,[(- 9 $ + c)iD^‘A>6x-^)] - ^Di,{-q<A> + c){D^mx-^) 

Airq Airq 

-^Dt{-q‘A> + c)x-^D^d<A> - ^Db{-q<A> + c)x-^D,mh^^ 

47rg Airq 


Vh ^ 

Airq 


+ ^(-g$ + c){x~^D‘^A>)SC^ 


bd 


= -^{D,<A>)iD^<A>)x-^Sx-^D\x-^D,<A>)6A> 

+ ^{Di,<^)x-^D,mh^^ - ^(4> - -)ix-^D^<A>)6C\d , (77) 

47r 47r q 


where we have used integration by parts twice and discarded the boundary 
terms. Since 1201 


1 


SC\d = -h^^DdShc, 


( 78 ) 
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we have 


(5L„ 


+ ^{Db<^)x-^D^^5h^^ 

47r 

- -){x-^D‘^^)Dd5hbc 
Stt q 

-^{D,^){D'>^)x-Hx-^D\u'^F,,)5^ 
47r 47r 

-^{D‘^^)x-^D^<^5hab 

47r 

'^,ab, 


+ —h^^Dd 
Stt 


($- -)x~ 
9 


(5/i, 


ab ■ 


(79) 


Hence, has been expressed as the linear combination of 6x^ and 
^hab- 
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